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TECHNICAL NOTE 1891 


ELASTIC BUCKLING OF A SIMPLY SUPPORTED PLATE 
UNDER A COMPRESSIVE STRESS THAT VARIES 
LINEARLY IN THE DIRECTION OF LOADING 
By Charles Libove, Saul Ferdman, and John J. Reusch 


SUMMARY 


Results are presented of calculations for the elastic Buckling load 
of a simply supported flat rectangular plate of uniform thickness subjected 
to unequal compressive stresses at two opposite edges, with a linear 
variation of stress between the two edges. The difference between the com- 
pressive stresses at the two loaded edges Is equilibrated by shear stresses 
along the other two edges. The results show that a plate with a linear stress 
gradient will buckle at an average stress that Is lower, but at a maximum 
stress that may be appreciably higher, than the uniform compressive 
buckling stress of the same plate. 

INTRODUCTION 


In most studies of the compressive buckling of flat rectangular plates 
the compressive stress has been taken as uniform throughout the plate. 

Cases of practical interest exist, however, in which the compression is not 
uniform but varies from one loaded edge to the other, as, for example, when 
the plate forms part of the upper or lower skin of an airplane wing in 
bending . 

One case of plate buckling under a load gradient is considered in the 
present paper. The problem is that of determining the elastic buckling 
load of a simply supported flat rectangular plate under a compressive 
stress that varies linearly from one loaded edge to the other, the 
difference in stress being equilibrated by uniform shears along the other 
two edges, (See fig. 1.) That this stress condition is a possible one, 
with regard to both equilibrium and continuity, is shown In the 
appendix. 

Solutions by the Rayleigh-Ritz method are obtained for several length— 
to— width ratios from 0 to », For each length— to-width ratio several 
ratios of minimum compressive stress to maximum compressive stress are 
considered. Some negative values are included among these stress ratios. 
These values correspond to a minimum stress that is tensile rather than 
compressive . 
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The influence of a loading gradient on the compressive buckling of 
a simply supported plate has been considered previously in reference 1. 

The loading gradient was exponential and the plate was assumed to have a 
thickness variation est ima ted to give minimum weigjit. Because economy 
of construction may often preclude the use of taper, the analysis of the 
constant— thickness plate in the present paper is felt to be of interest. 

The results of the present investigation are described in the section 
entitled "Results and Discussion." The theoretical derivations are included 
in the appendix. 


SYMBOLS 


a 

b 

P 

t 

E 

D 


plate length measured in compressed direction (x— direction) 
plate width 

length— to-^width ratio (a/b) 
plate thickness 
Young ' 3 modulus of elasticity 
Poisson's ratio 


plate flexural stiffness 


Bt3 

, 12(1 - ^)j 


y 


w 



coordinates of plate, shown in figure 1 
lateral deflection 

coefficients of Fourier series for w 


a 

0 max 


°min 


o 


av 


longitudinal compressive 

longitudinal compressive 
when buckling occurs 

longitudinal compressive 
when buckling occurs 

longitudinal compressive 
buckling occurs 


stress (stress in x— direction) 
stress at more highly compressed end 

stress at less highly compressed end 

stress at center of plate when 
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load— gradient parameter 


(5 


r mln \ 


max/ 


k = 


a tV 
« 2 D 


*%iax _ a ] 


max 


tt s 

it 2 D 




■^av ” a av 


tb g 

n 2 V 


dimensionless stress coefficients in terms of b 


t _ crta g 


k f = 


jt 2 D 


VI n 

max u max 


ta^ 

jt 2 D 


> dimensionless stress coefficients in terms of a 


k f = a 


ta c 


av av o 
it D 


potential energy 


°x> G y 


middle— plane direct stresses in x— direction and y— direction, 
respectively, prior to buckling, positive for tension 


*7 

i ) J, 

m, n, 

V, <1, r 


middle— plane shear stress in x-direction and y-^iirection prior 
to buckling 


> integers; also used as subscripts 


Primes ( f ) and double primes (") are used with subscripts i, m, n, 
and r to indicate odd and even values of the subscripts, respectively. 
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RESULTS AND DISCUSSION 
Graphical Presentation 


The results of the investigation are given graphically in figures 2 
to 5* In the first three of these figures the solid circles represent 
calculated points. The curves were obtained by fairing through the 
calculated points. There Is reason to believe that these curves should 
theoretically contain slight discontinuities In slope corresponding to 
sudden changes in the buckle pattern. The fairing of smooth curves, 
however, was felt to be Justified by practical considerations. 

The maximum compressive stress in the plate when buckling occurs is 
given in figure 2 for several different values of length— width ratio P 
where 

p-f t 1 ) 


and a is the plate length in the compressed direction and t>, plate 

width. The stress appears in the parameter k , defined as 

max 


■^max - a i 


tb 2 


max 0 
ji D 


( 2 ) 


where 

a 

max 

t 

D 

E 


compressive stress at the more highly compressed end when 
buckling occurs 

plate thickness 


plate flexural stiffness 


Et- 


12(1 - p 2 )y 


Young 1 s modulus 
Poisson f s ratio 


For each value of p, k^.^. Is plotted as a function of r, the ratio of 
the minimum compressive stress cT m ^ n to the maximum compressive 
stress crmax- This ratio Is usually a known constant independent of the 

absolute magnitude of the load. Negative values of r correspond to 
tensile values of cr min . 
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Figure 2 cannot contain a curve for p = 0 (b = 00 ) because the 
parameter k^^ involves b. The curves are therefore re plot ted in 
figure 3 in terms of k’^, where 

"Lj - 1 _ _ t a, / ^ \ 

max ~ u mx ”p~ l j ) 

nD 


Figures 2 and 3 show that the presence of a stress gradient (r ^ l) 
permits the plate to develop a that is greater than the uniform 

stress required for buckling (r = 1). For a given value of r this 
effect becomes less pronounced as p increases. The effect vanishes 
completely for P = ». 

Although the presence of a stress gradient permits the maximum stress 
in the plate to exceed its critical uniform value , it causes a reduction in 
the average stress at which buckling occurs. This reduction in the average 
stress is illustrated in figure 4- through use of the parameter k QTI ., defined 
as 


k 


av 


°av 



(«0 


where cTq,v the longitudinal compressive stress at the center of the 
plate when buckling occurs. 


The complete longitudinal direct-stress picture at buckling is shown 
in figure 5 for several values of p. The longitudinal compressive stress 
cr is plotted in terms of k for (3 other than zero and k T for p = 0, 
where 


gtb 2 

2 D 


(= 5 ) 


1 


k 1 


p 

cxta^ 

jt 2 D 


( 6 ) 


Figure 5 sho«r 3 at th© same tim© the reduction in the average critical 
stress and the increase in the maximum critical stress produced by a stress 
gradient, and how the increase in the maximum critical stress is diminished 
as the length— to— width ratio increases. 
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Tabular Presentation 

The basic calculated information from which the curves of figures 2 
to 5 were obtained is given in table 1. In the first part of the table , 
calculated values of k ay are tabulated against given values of r 

and 3 for the general case (3 ^ 0). The values of k mflT and k'^y 

plotted in figures 2 and 3 vere obtained from this basic information 
through the relationships 


^max 1 + r ^av 


k 'max - 


(7) 


( 8 ) 


In the second part of the table calculated values of k' av . for the 
infinitely wide plate (3=0) are tabulated against r, where k' ay is 
defined by the equation 


k' 


av ~ CT av 



( 9 ) 


The values of k’ ^y plotted in figure 3 for the infinitely wide plate 
were obtained from the tabulated values of k' aT through the relationship 


k' = —2 — k* 
max 1 + r av 


( 10 ) 


CONCLUSIONS 


A theoretical investigation has heen made of the buckling of a simply 
supported flat rectangular plate of uniform thickness under a compressive 
stress that varies linearly in the direction of loading. The parameters 
in terms of which the buckling stresses are plotted are the length—' to-width 
ratio £ and the ratio of minimum compressive stress to maximum compressive 
stress r. 
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The main conclusion to he drawn from the investigation is that when 
buckling occurs the stress at the more highly loaded end of the plate may 
he appreciably greater than the uniform buckling stress of the same plate. 

For a given value of r, this effect becomes more pronounced as P decreases. 
At the same time the average stress in the plate at buckling will be less 
than the uniform buckling stress. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., April U, 19^9 
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APPENDIX 

THEORETICAL DERIVATION 
General Case : 0 ^ 0 or 00 


E nergy expression .- Of all conceivable buckle patterns satisfying 
boundary conditions the actual buckle pattern is that for which the 
potential energy F, as given by the following expression, is a 
minimum (reference 2): 


F 


D 

2 





- 2(1 - u) 



S^w 

Sy 2 



>dx dy 





2r S$!Z 

dx 


5w 

dy 


dx dy 


(Al) 


where 

x, y coordinates of the plate , shown in figure 1 
w lateral deflection 

o x , Gy middle-plane direct stresses in x— direction and y-direction, 
respectively, prior to duckling, positive for tension 

t X y middle— plane shear stress in x— direction and y— direction prior to 

duckling 
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Middle— plane stresses .— Algebraic expressions for the middle— plane 
stresses cr x , g^., and T X y, needed in equation (Al), can be obtained by 

solving the plane— stress— equilibrium equations 


d g x + ^ T xy 


8x by 


= 0 


> 


C^CTy 

— £ + — = 0 

dx 


J 


(A2) 


First , the following expression for a x , corresponding to a linear variation 
of compressive stress between a mln and cr^^., is written 



(A3) 


where 


a 


av 


compressive stress in the x— -direction at plate center x 


_ g min = Minimum compressive stresg 
a max Maximum compressive stress 


a 

2 


a plate length 


Substitution of the expression for cr x from equation (A3) in equations (A2) 
permits these equations to be integrated for t xy and a y . This integration, 

subject to the conditions that no transverse compressive stress and no 
component of pure shear exist on the plate, gives the following expressions 
for T xy and a y : 




(AU) 


a 


y 


= 0 
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Through equations (A3) to (A5) the middle— plane stresses have all been 
expressed in terms of a single stress parameter a av and the load- 

gradient parameter r. 

Expressions (A3) to (Ap)* beside satisfying equilibrium conditions, 
also satisfy the condition of compatibility (reference 3) and, therefore, 
represent a physically possible set of stresses for an elastic sheet. 

Deflection function ,— The buckle pattern of the plate may be 
represented by the double Fourier series 


w 



N 

XL 


n=i 


w mn 


sln m* sin S2IZ 
a b 


(A6) 


which satisfies the boundary conditions for simple support term by term. 
Any arbitrary degree of accuracy Is In principle possible by letting M 
and N approach Infinity. For practical purposes, however, the use of a 
finite number of selected terms In equation (a6) Is sufficient. Terms 
that lie within the range m = 1 to M and n = 1 to N but that are 
not used may be thought of as having been assigned zero values for the 
coefficient v . 


The coefficients w^ are at present undetermined. Their relative 

magnitudes are subsequently chosen subject to the condition that F, as 
given by expression (Al), be a minimum. 


Stability determina nt.— Substituting equations (A3) to (A6) in 
expression (Al) and performing the indicated Integrations give 


r 


F « 


ff 2 P y < 
ab 


M N 


w, 


m=l n=l 


2 

mn "5“ 


£ (f + - 


nr 





mfcp = odd 


m^p 




2 


+ 2 


w w 
mn pn 
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where 



Ths deflection coefficienta w^ must now be chosen such that F is a 
minimum; that is, they must satisfy the equations 


SF 


Sw. 


= 0 


ij 


A = 1, 2, 3, . . 
\J = 1, 2, 3, . . 


. M 

. N j 


(A9) 


In expanded form, equations (A 9 ) can he written 


W *J T 


k av\P 


+ PJ 2 - i £ 


1 - r 


w , ig( l2 . + P 2 ) 
(i2 - P 2Y 
p±l=odd V ' 



1 + p=l PJ _ ^2 


h (^Afi YL w P q — ^3 

U + Pq (q2 _ j 2 )( l2 - p2) 

ttH = odd \ /\ ' 


= 0 


p+i = odd 
q±j = even 


(A10) 


i = 1, 2, 3, . . • M 

Sj = 1, 2 , 3 , • . . N, 


Equations (AlO) constitute a set of simultaneous equations to he solved for 
the w's. Since these equations are homogeneous, they can have nonzero 
solutions for the w's — that is, the plate can "buckle — only if the 
determinant formed hy the coefficients of the w's equals zero. This 
determinant contains kav- The lowest value of kav for which the 
determinant vanishes determines the buckling load. The order of the 
determinant is equal to the number of terms used in the deflection function 
(equation (A6)). 
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The accuracy of a solution thus obtained may be judged in two ways. 

The first consists in repeating the calculations several times, each time 
including additional terms in the deflection function (A6), and noting 
the sizes of successive changes in k ay . When these changes become 

negligible, the solution may be assumed to be converged. The correct 
buckle pattern, that is, the one for which F is a minimum, can also be 

shown to be the buckle pattern for which k^ is a minimum. Of several 

buckle patterns, therefore, the most accurate is the one giving the 

lowest k . The successive changes in k should therefore be negative, 
av av 

The second way of judging the accuracy consists in calculating the relative 
values of the coefficients w (their absolute values are indeterminate). 

The maximum w^ will occur for a particular set of values of m and n, 
say and n Q , and the magnitudes of the other will generally 

decrease as m and n become farther and farther from and n Q . 

The omission of any important terms from the deflection function (A 6) can 

be noted by systematically tabulating the relative values of the 
coefficients. 

Simplifications to stability determinan t.— The preceding paragraph 
describes in principle the use of the stability determinant to calculate 
buckling stresses. Two simplifications to the determinant are feasible 
and should be considered before calculations are made. The first 
simplification follows from the fact that the equations represented by 
equations (A10) fall into two completely independent subsets, one 
corresponding to even values of J and the other, to odd values of j. 

The first subset contains only those w f s for which the second subscript 
is even and the other subset contains only those w f s for which the second 
subscript is odd. The significance of the subscripts in the deflection 
function (equation (A 6)) indicates that the first subset corresponds to 

buckling that is anti symmetrical about the line y = — and the second, 

2 

to buckling that is symmetrical about that line. Calculations show that 
the second subset gives the lower buckling loads, and attention is herein- 
after confined to it. This subset is given as follows: 
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'iJ T 


if* *l*ftt{) £ 


P= 
p±i =odd 


ip(i 2 + p 2 ) 

( l2 _ p2^ 2 


- b 


M N 


1 ~ r \ y" w iJVg 

1 + r) p=1 q= f 3,5 M ( q 2 _ j2 )( i2 _ p2 ) 


= 0 (All) 


pti=odd 


'i = 1, 2, 3, • • • M> 

^ J = 1 , 3 , 5, • • • 


These equations. Instead of equations (A10), can now he used to write a 
stability determinant. 

A second simplification is possible by noting that equations (All) 
can be further divided into two subsets } one corresponding to even 
values of i and the other, to odd values of i. These two sets are 
written as follows > with primes on the odd subscripts and double primes 
on the even ones to make them more readily distinguishable. (Since J 
and q are always odd, the distinguishing mark is omitted from these 
subscripts.) . 


w i’j A i*j + 


t 


p"=2,4,6 


V’fi'p" 


M 


p"= 2 ,4,6 . . . q = 1,3,5 


N 

w „ C , „ =0 (A12a) 

f-5 p q i'Jp q 


( 


1* = 1 , 3, 5-, ... M 
<3 — 1* 5* • • * 
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and 


where 


w l"J A i"J + 


± 

=1 


w p'J B i ,, P* 


- £ 


N 

z 

. q-1,3,5 


w * C , if . 1 

p q i jp*i 


= 0 


(Al2b) 


i" = 2, 4, 6, . . . M\ 

, 3 » 1, 3, 5, . . . V 


A 


ij 




- 1 


2 




+ P 2 ) 

p 2 ) 2 


c ijp<i 



i.lpq 

(s 2 - J a )(i 2 -P^ 


Equations (A12a) and equations (A12b) are not Independent since the 
same w's appear in both sets; however, they afford the possibility of 
eliminating roughly half of the unknowns. For example, equations (A12b) 
may be solved for each coefficient of the type w m tt n in terms of all the 

coefficients of the type w m i n and the resulting expressions used to 
eliminate all the coefficients of the type w m ii n in equations (A12a). 

Equations (A12a) then become a set of equations involving only the coef- 
ficients of the type v « explicitly. The buckling stress is obtained 

as before by finding the lowest k av for which the stability determinant 
of these equations vanishes. If desired, ' equations (A12a) and (A12b) can 
be used to furnish an alternate set of equations in which only the coef- 
ficients of the type w « n appear explicitly. 
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The Just described, reduction of equations (A12a) and (A12b) to a 
single set of equations involving only the coefficients of the type w m , n 

is now shown in detail. Equations (A12b) are first solved for to 

obtain 


Wj II , = — 


i J 


M 

, Z 

P*=l,3,5, 


W P’J 


A i"j 


+ 


M 

n 

P'=l,3,5, 


N 


• q= 1,3,5, 

q^J 


W. 


p'q 


C l"jP'q 

A i"j 


(A13) 


/ i" = 2, k, 6, . . . M\ 
\ J = 1, 3, 5, • . . N/ 


In order to avoid confusion when equations (A13) are substituted in 
equations (A12a), the index q in equations (A13) is changed to r. 
Furthermore , in order to obtain a general formula for the coefficient of 
the type w m it n , the indices i" and J are changed to m" and n, 
respectively. The result is 


M 

z 

p =1,3,5, . 

Bm’V 

V , * 

p'n 

• • V’n* 



♦ ± 
P ,= l,3,5, . 

• 

ll 

• 

^W’np'r 
. . Wp ’ r V’n 

(A110 



• 

VO 

•s 

CM 

II 

. . m\ 



V n = 1, 3, 5, • 

. . N/ 
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In equations (A12a) the coefficients Vp«j and. w p"q are of the 

type v m t. n and may therefore he eliminated by the use of equations (Al4). 

Equations (A12a) then become 


M 

P -2, A, 6, 





VJ 


A p"J 


M W 

+ x x 

P ,= l*3,5, • • • **=1,3,5, 

rH 


v p f r 


C p"jP*r 

Vj 


X 

P" =2,4,6, . 



q.^J 





Bp" p « 


P q A "n 

P 1 


M 


P'=l,3,5, 



. r=l,3,5, 

r^q 


VpT r 


c Ap_'n 

Ap-’q 


= 0 


/i* = 1 , 3, 5, . . . m\ 
\ J = 1, 3, 5, • • . N/ 


or, by removing brackets 
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IT 


tf\ <r ’3 
* ^ 

on ct 1 


ft 

& 

c ft 

ft c 
o ft 

ft 


=WS' 


ITS 

on 

4C : 

• — l 

1 — 1 

1! 

II 

"ft 

"ft 


145 


a U > 
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A stability determinant is furnished by the coefficients of the v's in 
equations (A15). 

The main advantage of equations (A15) over equations (All) lies in 
the fact that, no matter how many w's are included in the deflection 
function (a 6 ), the order of the stability determinant will be only equal 
to the number of w's of the w m r n type (first subscript odd, second 

odd); whereas, if equations (All) are used they furnish a stability 
determinant the order of which is equal to the total number of w's. The 
main disadvantages of equations (AI 5 ) are that each term in the stability 
determinant is more complicated and that, if the relative magnitudes of 
t*he w's are desired, a separate calculation for w m it n coefficients 

(by using equations (Al4)) is required after the w m » n coefficients have 
been determined. 

Calculations .— Equations (All) rather than equations (AI 5 ) were chosen 
as the basis of a stability determinant. The matrix iteration method 
(reference 4-) was used to obtain the maximum value of l/k av (corresponding 

to the minimum value of and therefore the lowest buckling load) for 

which nonzero solutions for the w's are possible. The matrix iteration 
method at the same time furnished the relative magnitudes of the w's and 
thus permitted the judgment of whether or not all the important deflection 
coefficients had been included in the deflection function (A 6 ). The 
results of the calculations are summarized in table 2 which gives, for 
different assumed values of r and 3 , the calculated values of k ay and 

the relative magnitudes of the w^'b. 


Special Case: 3=0 

For the special case of an infinitely wide plate (b = « and 3=0), 
equation (A 6 ) has no meaning and the solution obtained is not valid. For 
this case, the plate will be assumed to buckle as a column, with straight 
lines in the y— direction. This buckle pattern can be represented by the 
single Fourier series 


M 

v = XI W m Bin 22 (Al 6 ) 

m=l a 


Equation (Al 6 ) is substituted in the energy expression (Al) in place of 
equation (a 6 ). The y integration in equation (Al) is now performed over 
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a unit width instead of from 0 to h. Minimization of the potential 
energy as before gives the following system of equations to replace 
equations (All) : 


1 8 Vk'av / Vl + r/ f=I 


vv 


p±.i -odd 


ip(l 2 + P 2 ) = 
(i 2 - P 2 ) S 


= 0 


(AIT) 


(i = 1, 2, 3, • . . m) 


The above system of equations was solved by matrix iteration to 
obtain k' av . and the relative magnitudes of the w m ' s used in 

equation (Al6). The results are given in table 3* 


Special Case; 3 = 00 

As for the case of 3 = 0, the buckle pattern of an infinitely long 
plate (a = 3 = 00 ) cannot be represented by equation (a 6). The solution 
for this case is, however, readily obtainable by physical reasoning. 

Consider a plate having a length a much greater than its width b, 
with a longitudinal compressive stress at "buckling varying from c m ^ n 

at one end to q at the other end. (See fig. 6.) The calculations 
for the general case (3^0) have shown that q mqT will be greater than, 
and q m1n less than, the uniform buckling stress of the plate, which in 

this case is . (See dash line in fig. 6.) The length of plate in 

b 2 t 

which the longitudinal compressive stress exceeds the uniform buckling 
stress is denoted by d in figure 6. 

Now assume that r = remains constant while a and, therefore, 

a max 

3 approach infinity. If remained constant, the length d would 

soon be many times the buckle wave length of the uniformly compressed 
plate. It is not physically plausible that a plate can sustain a stress 
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greater than its uniform "buckling stress over so great a length. Hence, 
for a given value of r, the following conclusion can be drawn! 


a max~^ 


t-rt 2 D 
b 2 t 


or 


^max - ^ ^ 


(A18) 


as p 00 . 
Also, since 


rk, 


■max 


+ k. 


■max 


av 


k 2(r + 1) 

av ' ' 


(A19) 


as P 


00 . 
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TABLE 1 


CALCULATED RESULTS FOR k nTr AND k' 

civ a v 


\ r 
p\ 

-1/3 

-1/5 

0 

!/5 

1/2 

1 

k 

a.v 

0.5 

4.55 

5.04 

5.58 

5.91 

6.16 

6.25 

.6 

3.66 

— 

4.53 

— 

5.05 

5.14 

.7 

3.13 

— 

3-9^ 

— 

4.44 

4.53 

1 

2.43 

2.77 

3.20 

3.53 

3.85 

4 


2.11 

— 

2.92 

— 

3.86 

4.49 

2 

1.93 

2.25 

2.71 

3.12 

3.63 

4 

3 

1.76 

2.07 

2.51 

2.92 

3.47 

4 

5 

1.63 

1.93 

2.36 

2.77 

3.33 

4 

00 

1.33 

1.60 

2 

2.40 

3 

4 

k av 

0 

0.823 

0.884 

0.941 

0.972 

0.993 

1 
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TABLE 2 — Continued 


CALCULATED VALUES OF v m AND CORRESPONDING VALUES OF k - Continued 

ELV 


3 

r 

■. 1 — 

n 




cm 






m = 1 

m = 2 

m = 3 

m = 4 

m = 5 

m = 6 





First Approximation 




1 

3 

1 

-0.00654 

-0.259 

-.OI56 

0.0297 

-0.00532 



3.21 


0 

Second Approximation 



1 

3 

5 



0.0296 

.00437 

-0.00518 



3.20 



First Approximation j 





-0.206 

-.0120 

0.0174 




3.54 

1 

1/5 

Second Approximation 





-0.209 

-.0120 

u 

0.0182 

.00265 

-O.OO336 



— 

3.53 



First Approximation 



1 

3 

1 

-0.123 

-.00672 

0.0061 

.00086 

-O.OOI5 



3.85 


1/2 




id Approximation 







O0U 01 






l 

3 

5 

1 

-0.0013 

-.00013 

-0.12 

-.0067 

0.006 




3.85 

V2 - 

-1/3 

l 

3 - 

1 

-0.0117 

— 0 . 607 
-.0102 

O.I76 

.0104 

-0.039 

-.0055 

mm 

m 


2.11 


0 

1 

3 ■ 



0.14 

.00724 

-0.024 C 
-.00343 

\OO58 

— 0 . 0016 

2.92 


^ACA^- 
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TABLE 2 — Continued 


CALCULATED VALUES OF v AND CORRESPONDING VALUES OF k - Continued 

XI1L1 aV 


P r 14 


m = l m = 2 m =3 m = 4 I m = 5 m = 6 m = 7 


1 1 -0.550 O.O7I5 -O.OO785 0.00201 -O.OOO518 

3 -0.00324 -.0040 .00313 -.00120 






-0.845 1 

.0102 -0.00103 


1 -O.798 1 

3 .OO887 


1 -O.718 
3 .OO70 


Jl -0.62 




1 -0.42 


1 O.380 
3 - .00469 




0.182 

.00430 

0.164 

.0066 


-0.48 0.135 

-.0067 .0050 


-0.44 0.11 

-.005 .0036 


-0.34 0.05 -0.01 

-.0035 .0016 


_j/s l 0.34 — 0.90 

' 3 -.004 .OO53 


1 0.27 -O.85 

3 -.0030 .0044 

^ 1 0 * 20 — 0 . 78 



3 -.0022 .0035-0.00028 -.0031 .0031 



1 0.11 -0.64 1 -0.50 

3 -.0011 .0021 -0.0002 -.0021 





-O.656 0.31 -0.118 

-.0048 .00571 

-0.64 0.29 -0.10 

-.0044 .0051 

-0.62 0.26 -0.084 

-.0037 .0041 


-0.59 0.21 -0.058 


-0.50 0.13 |-° .029 
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TABLE 2 — Concluded 


CALCULATED VALUES OF ^ AND CORRESPONDING VALUES OF k &y - Concluded 


p 


n 

v mn 


r 

m = 1 

m = 2 

m = 3 

m = 4- 

m - 5 

m = 6 

m = 7 

m = 8 

m = 9 

BL = 10 



-1/3 

1 

o.o8i*. 

-0.31 

0.66 

-0.97 

1 

-0.77 

0.48 

-0.26 

0.11 

-0.05 

1.63 


-1/5 

1 

0.07 

-0.27 

0.62 

-0.95 

1 

-0.76 

0.46 

-0.23 

0.1 

-0.04 

1.93 

5 

0 

1 

3 

0.053 

-0.22 

0.55 

-.0025 

-0.92 

.0021 

1 

-0.75 

-.0018 

0.43 

.0025 

-0,20 

0.072 


2.36 


1/5 

1 

3 

0.037 

-0.16 

0.48 

-.0019 

-0.89 

.0017 

1 

-0.73 

-.0015 

0.38 

.0020 

-0.16 

0.05 


2.77 


1/2 

1 

3 

0.018 

-.00026 

-0.083 

.00061 

0.33 

.0011 

-0.80 

.0011 

1 

-.00005 

-0.67 

-.0011 

0.28 

.0012 

-0.09 

-.00] 



3.33 
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Figure I- Plate subjected to unequal compressive stresses 
along two opposite edges. 
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Calculated points 
Faired curve 


Figure 3. Coefficient k / max used in calculating buckling 


stress <j. 


max* 


max 


= k' 

R maxVa 2 t ) 
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causing buckling , compared with uniform compressive 
stress required for buckling. 



